A Model for the Voltage Steps in the Breakdown of the Integer Quantum Hall Effect 
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In samples used to maintain the US resistance standard the breakdown of the dissipationless 
integer quantum Hall effect occurs as a series of dissipative voltage steps. A mechanism for this 
type of breakdown is proposed, based on the generation of magnet o-excitons when the quantum 
Ifall fluid flows past an ionised impurity above a critical velocity. The calculated generation rate 
gives a voltage step height in good agreement with measurements on both electron and hole gases. 
We also compare this model to a hydrodynamic description of breakdown. 

Pacs numbers: 73.43.-f 73.43. Cd 73.43. Lp 47.37.+q 
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In the integer quantum Hall effect (IQHE) regime [1], 
a two-dimensional electron fluid carries an almost dissi- 
pationless current and the ratio of the current, I^, to the 
Hall voltage, Vh, is quantized in units of e'^/h. How- 
ever, above a critical current, the dissipative voltage, Vx, 
measured along the direction of current flow, increases 
rapidly, leading to quantum Hall breakdown (QHBD). 
Several possible mechanisms for QHBD [2] have been pro- 
posed: avalanche heating [3]; percolation due to an in- 
crease of delocalized states [4] ; quasi-elastic inter-Landau 
level scattering [5-7]; acoustic phonon emission due to 
intra-Landau level scattering [8]; formation of compress- 
ible metallic filaments [9]; and resonant impurity scat- 
tering of electrons [10]. For certain samples, including 
those used to maintain the US resistance standard at the 
National Institute of Standards and Technology (NIST), 
breakdown occurs as a series of up to twenty steps, in Vx, 
of roughly equal height AI4 ~ 5mV, see Fig. 1 [11,12]. 
This is fundamentally different from the case in which 
breakdown is observed as a single large increase of Vx [3] . 

IQHE breakdown is not only of fundamental inter- 
est, but is also relevant to quantum metrology since a 
large value of Ix can improve the measurement precision. 
In this work we develop a theoretical model to account 
for the dissipative steps observed by the NIST group 
[11,12] and others [13,14]. We show that, in the presence 
of charged impurity-induced disorder, the quantum Hall 
fluid (QHF) is unstable when the local fluid velocity ex- 
ceeds a critical value. Under these conditions, magneto- 
exciton or electron-hole (e-h) pair excitations are gener- 
ated spontaneously near an impurity. The voltage step 
height, AVx, is directly related to the rate of formation 
of the pairs, which we calculate using a parameter-free 
model. This type of excitation of the QHF is analogous 
to vortex-antivortex pair formation in classical or quan- 
tum fluid flow around an obstacle [15]. 

Earlier work [16-18] has formulated a method for cal- 
culating the excitation dispersion relation of e-h pairs 
generated by exciting an electron from an occupied Lan- 
dau level n to an unoccupied Landau level (n+1), where 
n is the Landau level index. By extending these models 
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FIG. 1. Step like breakdown in the integer quantum Hall 
effect, from Ref. [12]. Each step contributes a multiple of 
~ 5mV to the dissipative longitudinal voltage. 

to include the electric field, arising from a charged im- 
purity and the Hall voltage drop across the sample, we 
obtain a critical electric field at which it costs no en- 
ergy to generate an e-h pair at a given wave- vector Q. 
We incorporate this into a calculation of the generation 
rate, W, of these pairs [19,20] due to a single charged im- 
purity. This then gives a dissipative voltage increment, 
AV; = hW/{2e). 

Our starting point is the Fermi golden rule, from which 
we calculate the rate of generation of e-h pairs, in the 
Landau gauge, due to a single charged impurity in a 
2DEG [19,20], 
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where V{r) is the impurity Coulomb potential, 
(j)n{r±, kx) is the electronic eigen- function in the xy plane 
and So(2:) is the envelope function of the first electronic 
sub-band. We assume that the n lower Landau levels are 
filled and the (n -I- 1) level is empty. There is thus no 
static screening of the impurity charge [21]. This gives 
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Here k is the dielectric constant (for GaAs k = 
O.llnFm"^). If we assume that the impurity is at the 
center of the subband then, 



M^n,(n+l) = y (^^2^ ) <5(en-e(„+l)) 
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where 
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^n{y — Yk^) is the simple harmonic oscillator solution 
to the Schrodinger equation centered on Y^^ = l%kx — 
Eym/{eB'^), m is the mass of the electron, Ey is the y 
component of the electric field, Ib = y/h/eB is the mag- 
netic length and B is the magnetic field. Calculating the 
transition rate out of state (n, k^) we find that 
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where [16-18] 
Ae = hu)c + 



f 1 k k' 

(6) 

-A„,(„+i)[A;^ - kx] - eEyl%{kx - kx) 
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and A„ [fc^ — fca;] includes the exchange and Coulomb 
local-field corrections, which arc independent of Ey. To 
obtain Eq. (7) we have assumed e^/naL ^ fi^jJc', hence 
we neglect Landau level mixing. In Eq. (6) there is an 
implicit condition for k'^ which must be calculated. This 
condition can be obtained from Eq. (7) when Ae = 0. If 
exchange and Coulomb interactions are omitted this con- 
dition is simply the same as for energy-conserving elastic 
inter-Landau level transitions. However, when interac- 
tion terms are included Q = k'^ — kx must be evaluated 
numerically. Eq. (7) can be split into two components: 
the excitation interaction energy, given by the first two 
terms, which is independent of Ey and the electrostatic 
energy, eEyl'^Q. In Fig. 2 the crossing point of the elec- 
trostatic energy (dashed line) and the excitation inter- 




FIG. 2. The magneto-exciton mode energy, 
fkiJc + /{4-kkIb)\)^ ,(i)„[Q] (solid line) for excitations from 



: to n = 1, for i5u 



0; the dashed line is the electrostatic 



energy, eEyl%Q. B = 12. ST and Ey = 1.5 x ICPV/m. 

action energy (solid line) gives the value of the e-h sepa- 
ration {1%Q) for which it costs no energy to generate e-h 
pairs, i.e. Ae = 0. To obtain the total rate of produc- 
tion of e-h pairs, we take Eq. (6) and sum over all initial 
states such that 



W 



2y^Wn = 2^ j dkxWn. 
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The above equation gives us the generation rate of e-h 
pairs by a single charged impurity at given B and Ey. 

Consider a local region of the sample where Ey is 
large enough to create e-h pairs, due to scattering from 
a charged impurity, at a rate given by Eq. (8). Such re- 
gions can be expected to occur at high current, possibly 
near the sample edge where the Hall field is expected to 
be large [22]. A pair created close to an impurity will 
drift along the Hall bar at a velocity v « EyjB, so one 
can imagine, for a fixed generation rate, a stream of e-h 
pairs moving along the Hall bar. Then, for j/ = 2, we 
have a situation in which a small fraction of electrons 
in the lower Landau level (n = 0) have been replaced 
by holes and the previously empty upper Landau level 
(n = 1) contains some electrons. As the c-h pairs move 
away from the high field region, the spacing between the 
electron and hole in a pair will increase and most pairs 
will eventually ionise by acoustic phonon emission. Due 
to the absence of empty states into which the excited elec- 
tron can relax and neglecting weak, second order Auger 
processes, we can assume that all the generated c-h pairs 
will eventually ionise and lead to a dissipative current 
% = eW, flowing across the Hall voltage equipotentials. 
At v = 2 this gives a dissipative voltage 



(9) 



We now compare the results of our model with break- 
down measurements which show voltage steps. The NIST 
experiments on the US resistance standard samples [12] 
were carried out a,t u = 2 and B = 12.37. The experi- 
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FIG. 3. The dissipative voltage, for electron [12] (a) and 
hole [14] (b) gas samples, calculated from Eq. (9) as a func- 
tion of Ey. The parameters used refer to the experimental 
conditions of [12,14]: (a) m* = O.OTme, = 3x lO^^m"^ and 
B = 12.3T; (b) m* = 0.4me, Ns = lO^^m"^ and B = 2.1T. 

ments show a series of dissipative steps in Vx of regular 
height AVx — 5mV, sec Fig. 1. In Fig. 3(a) wc plot our 
calculated dissipative voltage as a function of the back- 
ground electric field Ey. 

Since the rate of production of e-h pairs is strongly in- 
fluenced by the overlap between the wave- functions in the 
occupied (n = 0) and unoccupied (n = 1) Landau levels, 
it increases rapidly at a critical electric field. This critical 
electric field occurs when Eq. (9) gives a Vx compara- 
ble to the small background dissipative voltage governed 
by axx- The rate, and hence Vx, then reaches a max- 
imum when the electric field is such that the e-h pairs 
are formed close to the roton minimum of the magneto- 
exciton dispersion curve, QIb ~ 2.5. For a given Ey, 
for which the rate of production of c-h pairs is finite, 
the pairs formed at the breakdown point will relax, i.e. 
electrons in the upper Landau level relax their energy 
by moving towards one side of the Hall bar, whilst holes 
move in the opposite direction. This process tends to 
screen the Hall field over much of the Hall bar. Since 
the Hall voltage in the NIST experiments remains con- 
stant at its quantized value over the magnetic field range 
over which the dissipative steps occur, this screening ef- 
fect tends to enhance the electric field at the breakdown 
point. Thus, as the critical electric field is reached, the 
generation rate at the breakdown point increases rapidly, 
inducing a further increase in Ey, due to the screening of 
the Hall field in other regions of the Hall bar. For break- 
down at a single charged impurity we therefore expect 
the system to switch between two stable states, corre- 
sponding to 14 = and Vx ~ 5.6mV, which corresponds 
to the maximum value of Vx in Fig. 3(a). In the NIST 
experiments, a series of steps is observed, as seen in Fig. 
1 and we attribute each step to the formation of separate 
streams of e-h pairs generated by other charged impuri- 
ties, i.e. each step is associated with the local electric 
field at a particular impurity reaching its critical value. 
Fig. 3(b) shows the results of a similar calculation for 
the breakdown of the QHE observed in a hole gas sample 
[14]. In these experiments the step height was found to 
be AVx « ImV . Our calculation for this case gives a 
step height AVx = l-6mV, in qualitative agreement with 
experiment. 



The above calculated values are derived from a model 
which combines a calculation of the magneto-exciton dis- 
persion with an impurity-related tunneling rate. We have 
included explicitly the mechanism for the tunneling be- 
tween the Landau levels and exchange and Coulomb local 
field corrections. 

An earlier paper by one of us [15] drew an analogy 
between the process described here and the formation 
of vortices behind an obstacle moving relative to a fluid 
(e.g. the von Karman vortex street in classical hydro- 
dynamics). Using our model we now examine this anal- 
ogy more closely. One of the key ideas in understanding 
the quantum Hall states was the appreciation by Laugh- 
lin [23] that they have low Coulomb energy because the 
particles are tied to zeros in the many-body wavefun- 
tion. Zhang, Hansson and Kivelson [24] effected such a 
binding of zeros to particles using a Chern-Simons con- 
struction. This leads to an effective field theory for the 
fractional quantum Hall states in which the quanta of 
flux of a flctitious magnetic field are tied to each particle. 
This has two effects: flrstly, in the mean fleld approxi- 
mation the fictitious flux cancels the real flux through 
the system; secondly, the charge-flux composites become 
bosonic. Hence, at the mean fleld level, the QHF in a field 
corresponding to integer filling is replaced by a system of 
composite bosons in zero magnetic field. Such a system 
forms a charged superfluid and one can view the quan- 
tum Hall state as a composite Bose condensate. From 
this, Stone [25] formulated an effective superfluid hydro- 
dynamic model. The charged elementary excitations, e-h 
pairs, of the QHF appear naturally in this description as 
vortices in the order parameter for the composite boson 
superfluid. In this language, our model for the QHBD is 
the spontaneous creation of vortex-antivortex pairs when 
the QHF fluid velocity around an impurity reaches a crit- 
ical value [26]. 

According to Stone [25] the equation of motion of the 
QHF is 



m* [v - [v X Jl]] = e (E -h V X B) - V ( ^ |v|^ ■ 



(10) 



where v is the velocity field, /i is the local chemical po- 
tential containing all the interaction terms and CI is the 
fluid vorticity given by 

n, = n, = oandn, = 5^-^. (ii) 

ox oy 

Using Eq. (11), in conjunction with the continuity 
equation for the density of the QHF, and examining small 
pertubations in the velocity field, of the form 
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we find that when y/x = 0, 



Uw = Ae= —eEylBilsQ) + 



(14) 



This result is equivalent to Eq. (7) in the absence of 
interactions and for Ae = corresponds exactly to the 
elastic inter-Landau level tunneling condition [7]. Alter- 
natively within this hydrodynamic model it corresponds 
to the condition required to generate a vortex-antivortex 
pairs at zero energy. 

To make a direct comparison with our earlier quantum 
mechanical calculation we need to evaluate the dissipa- 
tive voltage drop along the sample due to the generation 
of these vortex-antivortex pairs from a single impurity. 
For a specific system we would have to rely on a nu- 
merical simulation of Eq. (10). However, we can make 
considerable progress by implementing what we already 
know about fluid mechanics [27]. Consider an obstacle 
in the path of a fluid. At a low fluid velocity the flow 
around the obstacle is laminar. When the flow rate is in- 
creased vortex-antivortex pairs are formed in the vicinity 
of the obstacle. However, a vortex street is not formed 
until the flow is fast enough to free the vortex-antivortex 
pairs from the local flow field near the obstacle. In this 
steady state each vortex-antivortex pair moves away from 
the obstacle at a velocity which is governed by the back- 
ground fluid velocity. This analogy suggests that the 
vortex-antivortex pair in a QHF moves away from the 
impurity at a velocity given by Ey/B. From classical 
hydrodynamics [27] it is also known that the distance 
between each vortex-antivortex (/) pair generated is ap- 
proximately three times the separation between a single 
vortex and antivortex (d) {d/l = 0.28). Now consider 
two states for our fluid: flrstly, the state where the fluid 
flow is laminar; then we know that the generation rate 
of vortex-antivortex pairs is zero. For the state where 
we have a street of vortex-antivortex pairs we have the 
condition Ae = 0, hence from Eq. (14) we can evaluate 
Ey/B. Dividing this velocity by the distance between 
each vortex-antivortex pair (d = Q.2SIb{IbQ)), to obtain 
the rate of generation of vortex-antivortex pairs, we can 
use this hydrodynamical model to estimate the voltage 
drop along the Hall bar as 

0.28^na>, 

To estimate the dissipative voltage step height from this 
hydrodynamic model we have to evaluate Eq. (15). At 
the moment we can only do this by referring back to the 
previous quantum mechanical calculation to give us an 
estimate of QIb- Taking the value for QIb for which the 
quantum calculations give the maximum value for Vx we 
find, using Eq. (15) that AT4; = 4.8/71^ for the electron 
gas {QIb ~ 1-9) and AT4 = 0.75my for the hole gas 
{QIb « 1.1). 

Above we have shown both a quantum calculation of 
the dissipative voltage drop along a Hall bar due to the 



presence of a charged impurity and an effective fluid 
model for QHF which also predicts the voltage drop along 
a Hall bar due to an impurity when a critical fluid ve- 
locity is reached. Both of these calculations predict volt- 
age steps which agree with experiment. We believe that 
the fluid model can be improved upon through numeri- 
cal modeling and the inclusion of interactions, however, 
the simple analysis provided above gives a strong clue 
that future work in this direction could be very useful in 
making direct comparisons with experiments. 
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